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Abstract 

The longitudinal dipole response of a quantum dot has been calculated in the 
far-infrared regime using local spin density functional theory. We have studied 
the coupling between the collective spin and density modes as a function of 
the magnetic field. We have found that the spin dipole mode and single 
particle excitations have a sizeable overlap, and that the magnetoplasmon 
modes can be excited by the dipole spin operator if the dot is spin polarized. 
The frequency of the dipole spin edge mode presents an oscillation which is 
clearly filling factor [v) related. We have found that the spin dipole mode 
is especially soft for even v values, becoming unstable for magnetic fields in 
the region \ < v < 2. Results for selected number of electrons and confining 
potentials are discussed. An analytical model which reproduces the main 
features of the microscopic spectra has been developed. 

PACS 73.20.Dx, 72.15.Rn 
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I. INTRODUCTION 



The far-infrared (FIR) response of quantum dots is a subject of current interest since the 
experiments carried out by Sikorski and Merkt |I| , and by Demel et al . These experiments 
and the subsequent theoretical work (see Refs. and references therein) have shown that 
the excitation spectrum of quantum dots in the FIR region is dominated by the dipole edge 
magnetoplasmon peak that splits into two different B dispersion branches when a magnetic 
field B is applied perpendicularly to the dot. These peaks are density (charge) collective 
modes excited by the operator Dp = Y^iLi^i- In the case of harmonic confinement by the 
potential |mu;Qr^, as a consequence of Kohn's theorem p the density mode is not coupled to 
any other mode, and the dipole operator Dp excites only two collective states at the energies 



uj± = \JuJq + jcOc i i'^c, where tOc is the cyclotron frequency. If the confining potential is 
not harmonic, Kohn's theorem does not hold, and on the one hand the energy of the modes 
depend on the number of electrons in the dot, and on the other hand a richer excitation 
spectrum appears. 

Raman spectroscopy has made possible to observe in the same sample single particle 
(sp), charge and spin density excitations |jlO|,|ll|, whose evolution as a function of B has 



been studied in recent experiments. This has revealed several interesting features of the sp 
|l^ and of the spin collective excitations [|13| in quantum dots. Limiting to the later, the 



experiments have determined that the spin mode lies very close in energy to the uncorrelated 
single electron excitations, and that magnetoplasmons can also be detected using spin- 
dependent probes. Besides, it has been experimentally determined that the spin mode has a 
much lower energy than the charge mode. These facts constitute the body of experimental 
results that any theory aiming at a quantitative simultaneous description of spin and charge 
density collective modes in quantum dots should reproduce. 

The dipole spin response function for unpolarized quantum dots at zero magnetic field 
has been recently addressed by two of us |T^. In the FIR regime, it has been found that the 
response is dominated by a low-energy collective dipole spin mode excited by the operator 
Dm = Z^ili Xi^l, where Xi and a* are Cartesian components of the position and spin vectors, 
and is the number of electrons in the dot. Similar modes have been described in atomic 



nuclei [^], and in alkali metal clusters |[T6|JTi 

The aim of the present work is to extend our previous study to the case of a quantum 
dot submitted to a perpendicular static magnetic field, which originates a B dependent spin 
polarization in the ground state of the dot. We will explicitly show that this not only causes 
the splitting of the spin dipole mode into two branches, one with negative and another with 
positive B dispersion, but also its coupling with the dipole density mode mainly excited 
by the operator Dp . We shall see that if the confinement is not harmonic and the dot 
is polarized, that operator also excites the dipole spin mode. Viceversa, when the dot is 
polarized, which is the case if it has an odd number of electrons, or for most cases when B 
acts on the dot, the spin response is coupled to the density response so that the external 
operator D^, = J2i ^i^^l also excites the density mode. When the system is fully polarized 
both modes coincide, while at zero polarization they are uncoupled. 

To this end, we have self-consistently evaluated the longitudinal response of the dot in 
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the time-dependent local-spin- density approximation (TDLSDA). As longitudinal we mean 
an external field which is either spin independent, or dependent on the spin component par- 
allel to the magnetic field, i.e., the z component. We present results obtained for selected 
numbers of electrons and confining potentials. Specifically, we have used a harmonic oscilla- 
tor potential to describe dots with = 5 and 25, and a disk confining potential to describe 
dots with 25 and 210 electrons, for which the FIR response has been determined in detail 
0]. The results obtained for the A^ = 5 dot have been presented as preliminary results in 
Ref. |T^. The ground state (gs) structure of the later two dots in intense magnetic fields 
has been recently addressed |]T^^. However, no self-consistent TDLSDA calculation for a 
dot as large as A^ = 210 has been carried out before even in the density channel. 

To obtain correct collective modes one needs to have a proper description of the ground 
state these excitations are built on. Several density functional calculations have addressed 
this question pT|-p5[]. The LSDA we have used in the present work is based on the exchange- 
correlation energy functional employed in Ref. PT[ as an input to construct their current 



density functional theory (CDFT). Within the range of magnetic fields we are interested in, 
we have checked that both LSDA and CDFT are yielding similar results for gs properties 
other than the current density. Tests of CDFT against exact and Hartree-Fock (HF) calcu- 
lations have been presented in Ref. pTI. Tests of unrestricted HF against exact gs energies 
in the filling factor region 2 > z/ > 1 have also been presented in Ref. |26] for small number 
of electrons (up to 5) . 

We conclude from the comparisons presented in the above references, that TDLSA may 
yield fairly accurate results for the density and spin response in the range of magnetic 
fields for which experimental information is currently available. Comparison with these 
experiments constitute the ultimate test of this essentially parameter-free approximation. 
Between 1/2 < v < 1, the results of Refs. P3| , p5[ seem to indicate that the exchange- 
correlation energy of Ref. |2l|] is not very suitable because it is based on the Levesque-Weiss- 
MacDonald interpolation formula , and improved exchange-correlation energy functionals 
like those of Ref. [^,0 are definitely better off (a detailed account of the accuracy of 
elaborated LSDA's in the < 1 regime is also given in Ref. |25|). It remains to be seen 
whether these improvements are crucial to describe the experimental data which for the 
moment are limited to comparatively small magnetic fields. 

To get a deeper understanding, simple methods have been developed to reproduce the 
gross features of the microscopic spectra. One such method is the sum rule approach used in 
Ref. []rU| to describe multipole density modes in quantum dots. We present here an analytical 
model, called vibrating potential model (VPM), which allows to understand the TDLSDA 
results in a clear way, thereby providing physical insight onto the longitudinal response of 
quantum dots. The VPM model has been widely used in nuclear physics p8[, where it was 
developed to describe nuclear collective modes. It has also been applied to the description 
of simple metal clusters p9|. 



II. THE LONGITUDINAL RESPONSE 

We consider a quantum dot made of A^ electrons moving in the z = plane where they 
are confined by the circular potential V~^{r) in the presence of a constant magnetic field B 
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in the z direction. In the local- spin- density approximation (LSDA), the single electron wave 
functions are given by the solution of the Kohn-Sham (KS) equations 



fair, 9) = ea(Pair,9) , (1) 



where = J df' p{f')/\f — f'\ is the Hartree potential. V^'^ = dSxc{p,rn)/dp\gs and 
y^xc _ Qg^^(^p^fy{^ I Qffi\^^^ are the variations of the exchange-correlation energy density 
Sxc{p,fn) in the local approximation taken at the ground state, and p{r) and m(r) are 
the electron and spin magnetization densities. The exchange-correlation energy density S^cc 



has been constructed from the results of Ref. on the nonpolarized and fully polarized 



two dimensional (2D) electron gas in the same way as in Ref. |2T|, i.e., using the von Barth 



and Hedin |^ prescription to interpolate between both regimes. 



We have used effective atomic units {h = e /e = m = 1), where e is the dielectric 
constant of the semiconductor and m is the electron effective mass. In units of the bare 
electron mass nif. one has m = m*me- In this system of units, the length unit is the effective 
Bohr radius = aoe/m*, and the energy unit is the effective Hartree H* = Hm* /e^. For 
GaAs we have taken e = 12.4, m* = 0.067, and g* = —0.44, which yields Og = 97.9 A and 
H* ~ 11.9 meV ~ 95.6 cm^^. In Eq. ([l|) Uc = eB/{mc) is the cyclotron frequency and 
Pb = eh/{2mf,c) is the Bohr magneton. The use of the same letter for the effective mass 
and the spin magnetization, and for the dielectric constant and the single electron energies 
should cause no confusion, since neither the mass nor the dielectric constant will explictly 
appear in the rest of the work. 

As a consequence of circular symmetry the cpaS are eigenstates of the orbital angular 
momentum i.e., (pa{r, 6) = M„^o-(r)e~*^^, with i = 0, ±1, ±2, . . .. The gs electron density is 
given by p(r) = J2a ri'a\ua{r)\'^ , while the gs spin magnetization density is expressed in terms 
of the spin of orbital a, {<Jz)a, as m(r) = J2ari'a{o'z)a\ua{r)\'^. The numerical calculations 
reported in the following have been performed at a small but finite temperature T < 0.1 K, 
and the KS equations have been solved by integration in r space. The thermal occupation 
probabilities Ua are determined by the normalization condition 

N = yn^ = y ^ (2) 

which fixes the chemical potential p. Our iterative method works for weak and strong 
magnetic fields as well, for which the effective potential is very different. It has probed to 
be very robust, allowing us to handle without any problem several hundreds of electrons. 

Once the gs has been obtained, we determine the induced densities originated by an 
external field employing linear-response theory. Following Refs. |]33| , |3^ , we can write the 
variation 5p„ induced in the spin density p„ (a =t, i) by an external spin-dependent field 
F, whose non-temporal dependence we denote as F = X)o- /cr(^) |o")(o"h 

5p„{r,u) = J2j dr'XcTa'ir,r';uj)f„,{f') , (3) 

a' 

where Xaa' is the spin-density correlation function. In this limit, the frequency 00 corresponds 
to the harmonic time dependence of the external field F and of the induced 5po-- Eq- ® is 
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a 2x2 matrix equation in the two- component Pauli space. In longitudinal response theory, 



F is diagonal in this space, and we write its diagonal components as a vector F = 
For the operators defined at the Introduction we then have 



fl 



X 
X 



and Dr 



X 
—X 



(4) 



The TDLSDA assumes that electrons respond as free particles to the perturbing effective 
field, which consists of the external plus the induced field arising from the changes produced 
by the perturbation in the gs mean field. This condition defines the TDLSDA correlation 
function Xo-o-' in terms of the free-particle spin-density correlation function x^Hl' through a 
Dyson-type integral equation: 



(5) 



cri(T2 



The free particle spin-correlation function at finite temperature is obtained from the KS sp 
wave functions, energies and occupation probabilities: 



rip 



-^<x,a'$]<^*(r)<^/3(r) . 

- + UJ + IT] 



(6) 



The label a (/?) refers to a sp level with spin a {a') and occupation probability Ua inp). To 
simplify the analysis of the results, we have added a small but finite imaginary part 7] to the 
energy uj. This will make an average of the strength function by transforming the 5-peaks 
into Lorentzians of width 27]. 

The kernel -ft"cro-' is the residual two-body interaction 



ri - r2 



+ 



d Pad Pa' 



6{ri - fs) 



(7) 



where 



d Pad Pa 



gs 



dp^ 



+ {r]a + r]a' 



gs 



dp dm 



gs 



dm? 



gs 



= K{r) + {r]a + i]a') L{r) + rjaTja' I{r) 



(8) 



with 7]^ = 1,?7| = —1. The last expression is the definition of the K, L, and / functions. 

When the system is not polarized, there are only two independent correlation functions. 
These are Xpp ^-nd Xmm describing, respectively, the density response to Dp and the spin 
response to Dm- They are given by 



Xpp — Xtt + Xll + Xu + XiT 

Xmm = Xn + Xll ~ Xtl ~ Xll 



(9) 



and the four equations (||) reduce to two uncoupled equations for Xpp and Xmm whose kernels 
are given by l/ri2 + K6{ri2) and /5(ri2), respectively, and the free-particle correlation 
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function x^"^ = x^i' + Xu = "^X^i same in both channels because Xp;' = Xu- This 

constitutes the paramagnetic hmit of the longitudinal response with uncoupled density and 



spin channels |]T4[, in which the residual interaction consists of a Coulomb direct plus an 
exchange- correlation terms in one case, and only of an exchange-correlation term in the 
other. 



(0) 



When the system is polarized one no longer has x\i' 
independent correlation functions 



(0) 



and there are two more 



Xpm = xtt ~ Xii - Xn + XiT 

Xmp = Xn - xii + xti - xu , (10) 

which produce the density response to Dm and the spin response to Dp, respectively. 

Equations (^) have been solved as a generalized matrix equation in coordinate space 
after performing an angular decomposition of Xaa' and K^„' of the kind 



9') 



(11) 



Only modes with 



±1 couple to the external dipole fields Dp and Dm- This can be 



readily seen performing the angular integral in Eq. (|^). In practice, we have considered the 
multipole expansion of the external field, using the dipole vectors 



-re 



and 



-re 



(12) 



For a polarized system having a non zero magnetization in the gs, the i = ±1 modes are 
not degenerate and give rise to two excitation branches with AL^ = ±1, where is the gs 
orbital angular momentum. 

The response functions corresponding to the above dipole fields have been obtained from 
the £ = ±1 components of the correlation functions Xab (^) ''"'j ^) with A, B = p,m as: 



MB V^; "T ^AB (13) 

Their imaginary parts are related to the strength functions as Sab{!jj) = ^^^[c^ABi^)]- 
Although the excitation energy u and strength functions are always positive, it may be 
easily verifed that the following relations formally hold: 



Re 
Im 



O^AB l^J 



Re 



a 



-UJ 



a 



{-£] 

AB 



— Im 



AB\ 



(14) 



To check the numerical accuracy of the calculations we have used the f-sum rules for the 



dipole operators, which can be expressed in terms of gs quantities |35 



Spp{uj) ujduj 



-(0|[D„ [H,Dp 



N 
~2 



m 



{mm) 



m 



(mp) 



(pm) 

m{ 



1 

{uj)ujduj = -(0|[D^, [H,DmW) = - 



j Smpit^) ujduj + J Spm{i^) ujduj = (0| [H, Dp]]\0) = 25*2 , (15) 



where Sz is the total spin of the ground state. 
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III. RESULTS 



Figures display the dipole strength function of the = 5, 25, and 210 dot for 
selected B values. Solid lines correspond to the density response to Dp, and dotted lines to 
the spin response to D^, that is to Spp and to Smm- Dashed lines represent the free particle 
strength function. For the five electron dot we have used a parabolic confining potential 
V~^{r) = ^mulr"^ with uq = 4.28 meV, and for the other dots we have used a disk confining 
potential as described in Refs. [0,0 . 

For = 25 and 210, most B values displayed correspond to integer filling factors u. It 
was found in Ref. |^ that for = 210 and an R disk confining potential, these values follow 
the law u = 27ccns/{eB) pertaining to the 2D system, where Ug = N/{tiB?) is the electron 
surface density. For the = 25 dot, that law yields B = 3.29 T as the value at which the 
system becomes fully polarized. Actually, we have found that the dot is in the maximum 
density drop (MDD) state for 3.42 T < B < 3.70 T. We have taken as u = 1 configuration 
that corresponding to Bi = 3.56 T, and have defined the other u configurations as those 
corresponding to the value B^, = Bi/u. 

The results of Ref. for the N = 25 dot seem to indicate that a confining potential of 
parabolic type might be more adequate (see the discussion of the charge mode at the end 
of this Section). Consequently, we have also studied it using a harmonic confining potential 
with ujq = 2.78 meV to reproduce the experimental dipole energy at i? = 0. In this case the 
system is in the MDD state for 4.46 T < B < 4.69 T. We have taken Bi = 4.58 T and have 
defined the other u configurations as indicated before. 



The existence of a MDD state for this dot is in agreement with the findings of Ref. p4 
It is also worth to remark that even for such a small dot, u as defined before also coincides 
with the number of occupied (n, T or |) bands for values up to z/ = 5-6, which correspond 
to rather low B intensities. We present in Fig. ^ the strength function corresponding to A^ 
= 25 with parabolic confinement. 

Figures |1]-^ show that in both channels the response at -B = is concentrated within 
a small energy range, with one single peak or with several closely lying fragmented peaks 
which exhaust most of the f-sum rule. The peak energy is lower in the spin than in the 
density channel. This is due to the character of the residual interaction, which is attractive 
in the former and repulsive in the later channel, and shifts the TDLSDA responses from the 
free particle response in opposite directions. The residual interaction in the spin channel is 
weaker than in the density channel, where not only the exchange-correlation but also the 
Coulomb direct term contributes. Consequently, the spin response is close in energy to the 
free response. It is thus difficult to distinguish the collective spin mode from the single 
particle spectrum. In large dots it also causes a stronger Landau damping in the spin than 
in the density channel. These facts have been observed and discussed in Ref. [p!3| . 

At i? = 0, as a consequence of Kohn's theorem, if the confining potential is harmonic 
with frequency Uq the excitation energy of the dipole density mode is equal to Uq irrespective 
of the number of electrons. Otherwise, the excitation energy depends on A^, see for example 
Refs. P,p!4|. In the spin channel Kohn's theorem does not hold, and a size dependence 



appears in the dipole spin mode even for parabolic confining potentials |jT4|. 

When B is not zero, the dipole mode in either channel splits into two branches, one with 
negative and another with positive B dispersion. The splitting is due to the breaking of the 



7 



^-degeneracy of the sp energies by the apphed magnetic field. Several new phenomena then 
appear. We first notice that for B values such that the spin of the dot gs is different from 
zero, the spin and density modes are coupled. This is particularly apparent in the = 210 
dot. Indeed, at -B = 1.71 and 5.14 T the system is almost paramagnetic, having 25*2 = 2 
and 0, respectively (see Fig. 4 of Ref. [^). As a consequence, the modes are uncoupled, as 
it can be seen from panel (b) and especially from panel (d) of Fig. |^. In contradistinction, 
ai B = 3.43 and 7 T we have 2Sz = 54 and 74: the system has a large spin magnetization in 
the gs and the spin and density modes are clearly coupled, as displayed in panels (c) and (e) 
of that figure. One can see a distinct peak in the spin response at the energy of the density 
mode. This effect has been experimentaly observed . The strength of this peak increases 
with Sz and when the system is fully polarized, which happens sligthly above B = 10 T 
for the iV = 210 dot, all the strength is transferred from the spin to the density channel. 
Conversely, the spin mode may be observed in the density channel with some intensity. This 
effect is hindered because Kohn's theorem prevents it from occuring for parabolic potentials, 
and for the disk potential it is of order {23^ /NY. 

The interplay between charge and spin modes is especially marked in the mixed channel 
where the density response to the spin dipole operator Dm, or the spin response to the 
density dipole operator Dp are described. This is shown in Fig. ^ for the N=210 dot ai v = 
3. One clearly observes two peaks at the energy of the density modes, and another two at 
the energy of the spin modes. It can be understood casting the mixed response into a sum 
over 'spin dipole states' |m) and another over 'charge dipole states' \p) 

Smp{uj) = Spm{uj) = '^{0\Dp\n){n\Dm\0)6{u - Uno) 

n 

= J2{0\Dp\p){p\Dm\0)6{uj-Upo) 
p 

+ J2{0\Dp\m){m\Dmm{uj - iOmo) • (16) 

m 

For a disk confining potential, the matrix element {0\Dp\m) is not zero and there is a 
contribution to Smp from the spin modes. For a harmonic confining potential {0\Dp\m) is 
zero and only the density modes would contribute to Smp through the p-sum in Eq. (|16|). 

The B dispersion of the main peaks of the spectrum is reported in Figs. |^— In these 
figures the density modes are represented by dots and the spin modes by triangles. Solid 
symbols correspond to integer filling factor values, and the inserts show the negative B 
dispersion branch of the spin mode. As a guide, we have drawn lines starting at the value 
of the B = frequency and following the VPM B dispersion laws (see Eqs. 

Several features of these figures are worth to discuss. Concerning the spin modes, we 
first see that at low B their energy is much smaller than the energy of the density modes, 
in agreement with the experimental findings of Ref. [|1^]. At higher B the dot is eventually 
fully polarized and the longitudinal spin and density modes merge, as in the two dimensional 



electron gas (see Fig. 5 of Ref. ||36|). This is not explicitly shown in the figures. Secondly, 
the negative B dispersion branch of the spin mode manifests a clear oscillatory behaviour 
with z/, similar to that found for the density response p^, also discussed in Ref. [BO]: the 



'paramagnetic' even u configurations have softer spin modes than the 'ferromagnetic' odd u 
configurations. 
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Our calculation predicts a spin instability occurring when the energy of the spin mode 
lower branch goes to zero at a critical B between u = 1 and 2. This instability also manifests 
in the static spin polarizability Re[a^^''(0)], which becomes negative at these large S's. This 
indicates that the gs is no longer an energy minumum and thus is not stable. It is worth 
to recall that no collective spin dipole modes are observed in the experiments at these 
rather high B values. However, one cannot discard this might be due to the strong Landau 
damping existing in this energy region. We remark that within LSDA, a spin wave instability 
of similar type has been found in parabolically confined dots at i? = |^ whose origin is 
a spontaneous breaking of the circular symmetry, which could also be the origin of the spin 
dipole instability. 

It is of some interest to look at the effect that a moderate temperature may have on the 
spin response of the = 210 dot at high magnetic fields. As an example, we have considered 
5 = 5.14 T (i/ = 2), and B = 7 T, and have obtained the spin and charge responses at T 
= 2 K. 

The results are displayed in Fig. 0. Apart from the well known T-independence of the 
charge response at low temperatures, the figure shows that the low energy spin peak recovers 
at T = 2 K, no longer lying at zero energy [Q. We do not display the results for u = 2, 
since the response remais essentially unchanged. The different behavior is due to the quite 
distinct sp spectrum of these configurations, as can be seen from Figs. 2(b) and 5(d) of Ref. 
||20| , which makes the B = 7 T state very sensitive to thermal effects (see Fig. 7 of Ref. 
20|). It is a general rule that integer u states are rather inert thermally because for them 



the chemical potential lies between Landau bands and consequently, very few sp states are 
affected by a change in T if it is much smaller than the sp energy gap, which is the present 
case. The situation is completely different for B values that do not correspond to integer u 
values. In this situation, the chemical potential is on top of a Landau band, and the changes 
are sizeable because of the large density of states around the chemical potential (see Figs. 
2, 5 and 7 of Ref. |2g). 



Finally, we would like to comment on the density dipole mode. For the parabolic confining 
potential. Figs. ^ and |^ show the well-known result that the density response yields the 
classical law represented by the first of Eqs. |2^. For the disk confining potential, that law is 
also fairly obeyed, particularly by the negative B dispersion branch. We have systematically 
found that the positive B dispersion branch is fragmented, especially for the A^ = 210 dot. 
Comparing our calculations for the N = 25 dot with the experimental results , we conclude 
that a parabolic potential is better suited than a disk potential to represent the physical 
situation. On the contrary, the confining potential of the A^ = 210 dot is not parabolic, and 
this is the origin of the second upper branch found in the experiment. The non-harmonicity 
of the confining potential has been presented in Ref. as the origin of that branch on the 
basis of a Hartree+random phase approximation (RPA) calculation in dots with A^ < 30 
(see also Ref. |^). Our TDLSDA calculation supports that interpretation. 

It is worth to see that the density response of such a large dot displayes two instead of 
one satellite branch (compare with the results for the small N = 25 dot), i.e., TDLSDA 
seems to produce a high frequency peak that is more fragmented than in the experiments. 
We want to remark that the upper branches disappear at intense B, in agreement with the 
experimental findings. This gives further support to our explanation and that of about 
the origin of these branches. Indeed, at higher magnetic fields one expects that the harmonic 
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cyclotron potential dominates over the other contributions in the KS equation. 



IV. THE VIBRATING POTENTIAL MODEL 

The intuitive idea behind TDLSDA is that a small amplitude time-dependent variation 
of the mean field around the static equilibrium configuration produces an oscillation in the 
electronic density, which causes a small amplitude collective motion of the system. This 
motion is self sustained if the induced density is precisely that needed to generate the 
oscillating potential. The vibrating potential model naturally arises when one considers the 
first iteration of the (perturbed mean field) < — > (induced density) selfconsistent scheme. In 
homogenous systems, where translational invariance determines the shape of the induced 
density, it yields the exact solution. In finite systems (nuclei, metal clusters, dots), the 
model provides a useful approximation. 

Using the general method described in the Appendix, we consider the following VPM 
Hamiltonian 



N 

i=l L 

N 

= ^ [Ho{f,) + SVir,,t)] 
1=1 

with a time-oscillating potential 



(17) 



5V{r,t) = -- 



^Xk)x - 2 — 

k 



k 



(18) 



where Uq = nng/R. The time dependence of 6V{f,t) is in the (...) spatial foldings with 
the densities induced by a time dependent external field, see the Appendix. To get the 
static part Hq of the Hamiltonian Eq. (|T^ we have assumed an exact cancellation between 
the Hartree and the external potentials, and have taken a parabolic approximation for the 
exchange- correlation potential at i? = 0. is the exchange-correlation constant introduced 
in the Appendix. 

This VPM Hamiltonian can now be solved analytically within RPA by finding the oper- 
ators solution of the equations of motion 



(19) 



We have used the methods illustrated in Ref. 
Hamiltonian as 



to compute the conmutators with the 



[H,0] = {H,,0] + 5V{0) 



(20) 



where Hq is the static Hamiltonian and 5V{0) the variation arising from the induced den- 
sitites. It can be shown that the solutions to Eq. (pj]) are given by 
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N 



N[l 



(21) 



where 



N N 

Q± = =t Vi) ' P± = J2iPxi ± Pyt) 
i=l i=l 

N N 



(22) 



i=l 



1=1 



and 



\ 



(23) 



The corresponding frequencies are 



and it is easy to verify that 



±0 



± 



(24) 



(25) 



The states \uj±'^) = 0±"^ |0) are orthonormal and carry an orbital angular momentum 
Lq ± 1 and a spin 5*0, where Lq and 5*0 are the orbital and spin angular momenta of the 
ground state, respectively. 

The charge dipole and spin dipole strengths are distributed among the above states as 
follows: 



P\|2 



P\|2 



P\\2 



IN 
4 to 

C2 



D,|O|^ = |(0|D,|c.")|^ = 



m\ |2 



rn\ |2 



1 N 



1 



V N 



21 



(26) 



It is a simple matter to check that the above solutions exhaust the sum rules Eq. (pTSl). 

This vibrating potential model reproduces the most salient features of the full self- 
consistent calculation. Its parabolic form guarantees that Kohn's theorem is fulfilled, and as 
a consequence the response to Dp is shared by just two peaks, which have the same strenght 
N/Auj. Accordingly, the spin dipole modes leu™) are not excited by the dipole operator Dp 
and the corresponding matrix elements vanish. Another consequence of Kohn's theorem is 
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that within VPM only the density mode contributes to the mixed response. The model 
predicts that in the mixed channel, the charge dipole modes \uj±) are excited with the same 
strenght Sl/uN by the spin dipole operator D^. 

In the spin channel, the spin dipole operator excite both the charge \uj^) and 
spin modes. These peaks have strengths | (OlD^lu;^) p = A^[l - {2SjNf]/Au)„ and 
I (0|Z)m|i^± ) P = S1/ujN , which are the same for both AL^ = ±1 branches. Finally, when 
5*2 = the density and spin modes are decoupled, and when the system is fully polarized, 
i.e., 2Sz = N, all the strength is transferred to the density channel. 

Besides this qualitative agreement, the strengths given by the first and last lines of Eq. 
(p6|) agree well with the result of the microscopic calculation. Also the ratio {2Sz/NY 
between the strength of the density peaks excited by Dm and Dp is reproduced. 

The second Eq. ( p^ can be used to determine the energy of the spin dipole mode a;™ 



as the first one has been often used in the case of the charge dipole mode, if we fix the 
parameters entering that equation. We first take uj^^ = 27ms/ R"^- This estimate is obtained 
identifying the kinetic energy per particle in the exchange-correlation harmonic oscillator 
potential with that of the 2D Fermi gas 7ms/2, and approximating the mean square radius 
of the dot (r^) by R'^/2 = rfA^/2, where = 1/(7™^). One then gets 



irW. + b + 4±^ ■ (27) 



\ rfiV' « " ' 2' 4 2 



The value of r^V^ is related to the spin susceptibility of the two-dimensional electron gas 
Xo/x = rl K + 1 (Ref. H). Eq. (|7D yields values for uo"^ which agree with the TDLSDA 
ones. In particular, we have checked that the dependence of the 5 = calculations 
reported in Ref. is well reproduced taking r^\4- — —0.3 at ~ 1. 



Notice that the A^ and dependence of energy of the spin mode is different from that 
of the density mode, which is given by 



' (28) 



\ r3ArV2 4 2 



when we take for the frequency ujq of the confining potential the estimate c^o = r^^/^A^"^/^ 
obtained from the disk potential in the r/R « 1 limit. 

We report in Fig. |TT] the energies given by Eq. (P7|) for a dot of A^ = 200 together with 
the experimental spin dipole mode of Ref. |T^. A value of = 0.65 has been used which 



yields r^V^ = —0.24. From the figure one sees that the negative B dispersion experimental 
branch is fairly well reproduced. It has been experimentaly found that the dipole positive 
B dispersion branch, as well as other spin modes laying at higher energies are rather B 
independent, so one should not expect that the simple model leading to Eq. (^) reproduces 
this behavior. 

The VPM works better for large than for small dots. This is not surprising in view of 
the approximations leading to VPM. It can be seen comparing with the TDLSDA results 
in Figs. after being advised that for the A^ = 5 dot, the dotted line representing \uj'^) 
is not passing near the main peak energies of the positive B dispersion branch, but near to 
the satellite peak energies. 
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Eq. p7| ) shows that within VPM, the spin instabihty at i? = mentioned in the previous 
Section occurs when r^V^ + l/2<0. In the 2D electron gas it happens when + 1 < 
30| . The difference between these conditions is due to finite size effects, which are crucial 



to determine the values of at which the instability appears in quantum dots. Whereas 
in the bulk the spin instability sets in at ~ 37, Eq. ( p7D yields that in dots it happens 



at Ts ~ 3. This value is well within the range of those found in Ref. |3T[ for the onset of a 
spin density wave instability in small magic dots, and a factor of two larger than the typical 
values obtained in that reference for open shell dots. 

V. SUMMARY AND OUTLOOK 

In this work we have thoroughly discussed the longitudinal dipole response of quantum 
dots. We have shown that TDLSDA is able to reproduce the main features of the exper- 
imental results. In particular, we have found that the density and spin modes are clearly 
coupled in the spin channel if the system is partially polarized, and that the spin modes are 
especially soft for even filling factors. We have determined that at 5 7^ the spin response 
lies in the energy region of single electron excitations. 

We predict that the frequency of the spin edge magnetoplasmon presents an oscillatory 
behavior as a function of z/, and that the spin dipole mode becomes unstable at B values in 
the 1 < u < 2 region. 

Our numerical scheme has allowed us to study large size dots whose spectrum has been 
experimentally studied in detail, instead of relying on extrapolation of the results obtained 
for small size dots. This is crucial to identify the u behavior of physical quantities like 
excitation frequencies. 

TDLSDA can be easily applied to other multipole spin and density excitations. This is 
relevant in view that recent experiments have been able to identify the spin monopole and 
quadrupole modes, and likely charge modes different form dipole |jl3l. Work to extend the 



present study to other multipolarities, and to study the spin transverse channel in large dots 
along the line described in Ref. [0 is under way. 
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APPENDIX: 

In this Appendix we provide a general derivation of the VPM Hamiltonian leading to 
Eq. (p!?!). We start from the density and magnetization variations 

{Sp{f,t),6m{r,t)) = a(t)(y xPo{r) ,V xmo{r)) 

+ /3(t)(V,mo(f),V,po(r)) , (Al) 
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obtained from the static ground state density po and spin magnetization mo through the 
unitary transformations e"^*'' and e'^*-*^ '^I'^l ^ respectively, in the hmit of small defor- 
mations. In Eq. ( |A1| ) a{t) and (3{t) give the amplitude of the oscillations. 

Assuming for the exchange-correlation potentials the forms V^'^ = V^'^{p) and W^"^ = 
mVcr/po, which amounts to expand S{p,m) around m = up to order and identifying 
Va with pod'^Sxc{po,iTi)/d'm'^\m=o, one obtains the variation in the one body potential of Eq. 
(|I|) induced by the density variations (|A1|): 



25 

6V{ri, t) = a{t) [ Q(rl) + ^Q^lf.) al ) + P{t) 



2^ 

A^ 



Q{ri) + Qair^) crl 



where 



(A2) 



Q{r[ 



■dr' 



Qa{r) = — V^po(^" 
Po 



(A3) 



with Vq = V^'^{p = Po)- We have further assumed that mo = ^Po in the gs. Using the 



results 



{J2Q{f^) ^ J Q{f)6p{f,t)df= - (a(t) + '^m) J Po{r)VMr)df 



(EQ<x(rl)a;) = J Q,{r)5m{r,t)dr 



25, 
A^ 



a{t) + I3{t) / po{r)V.Q.{r)dr (A4) 



it is then possible to write the variations in the one-body potential of Eq. (0) in the separable 
form 



25. 



(E.Q(r1)) 



25, {Y.iQa{r^)(yl) 



J(V..g(f))po(f)rff N J{VMr))po{r)df^ 
' (E.g.(rl)a:) _ 25, (E^Qlf,)) ^ 
J(V.Q.(f))po(f)df N J{VMr))Po{f)dr, 

'25 



(A5) 



Following Ref. [^], one could now express the various responses to an oscillating potential 
of the form ( [A5D in terms of the independent particle response functions through RPA- 
type equations. However, our aim here is to develop an analytical model which allows us to 
understand in a simple way the numerical results of Section III. To this end we take in (|A3|) 
a harmonic oscillator for the one-particle potential Vq = to simulate the short range 

effects, and a step function for the electronic density po entering the long range contribution 
/ Va;Po(^)/|^— r'\dr'. We also assume Qaifi) ^ —^VaXi/ (r^). Equation ( p!?] ) of Section IV 
is then obtained. 
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FIGURES 



FIG. 1. Dipole strength function (effective atomic units) of the N = 5 dot as a function of 
frequency (meV). Sohd and dotted hnes correspond to the density response to Dp, and to the spin 
response to Dm, respectively. Dashed hnes represent the free particle strength function. 

FIG. 2. Same as Fig. 1 for = 25. 

FIG. 3. Same as Fig. 1 for AT = 210. 

FIG. 4. Same as Fig. 2 using a parabolic potential with loq = 2.78 meV instead of a disk 
confining potential. 

FIG. 5. Mixed Smp{oj) response function (effective atomic units) of the N = 210 dot at = 3. 
The arrows indicate the density and spin mode peaks. 

FIG. 6. B dispersion of the main peaks of the N = 5 spectrum. The circles correspond to 
density modes and the triangles to spin modes. The solid symbols correspond to u = 1. The lines 
represent the VPM B dispersion laws with fitted value at B = 0. 

FIG. 7. B dispersion of the main peaks of the = 25 spectrum for a disk confining potential. 
The circles correspond to density modes and the triangles to spin modes. The crosses are experi- 
mental points from Ref.[2]. The lines represent the VPM B dispersion laws with fitted values at 
B = 0. The insert shows the negative B dispersion branch of the spin mode. From left to right 
the solid symbols correspond to = 6 to 1. 

FIG. 8. B dispersion of the main peaks of the A^ = 210 spectrum for a disk confining potential. 

The circles correspond to density modes and the triangles to spin modes. The crosses are experi- 
mental points from Ref.[2]. The lines represent the VPM B dispersion laws with fitted values at 
B = 0. The insert shows the negative B dispersion branch of the spin mode. From left to right 
the solid symbols correspond to = 8 to 1. 

FIG. 9. Same as Fig. 7 using a parabolic confining potential with uo = 2.78 meV instead of a 
disk confining potential. 

FIG. 10. Dipole responses for the A^ = 210 dot at S = 7 T. Top panel, T = 0.1 K. Bottom 
panel, T = 2 K. Lines as indicated in Fig. 1. 

FIG. 11. B dispersion of the spin dipole mode within VPM. Experimental points are from Ref. 
[13]. 
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